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Abstract 

We have employed non-adiabatic molecular dynamics based on time-dependent 
density-functional theory to characterize the scattering behaviour of a pro- 
ton with the Li4 cluster. This technique assumes a classical approximation 
for the nuclei, effectively coupled to the quantum electronic system. This 
time-dependent theoretical framework accounts, by construction, for possi- 
ble charge transfer and ionization processes, as well as electronic excitations, 
which may play a role in the non-adiabatic regime. We have varied the in- 
cidence angles in order to analyze the possible reaction patterns. The initial 
proton kinetic energy of 10 eV is sufficiently high to induce non-adiabatic 
effects. For all the incidence angles considered the proton is scattered away, 
except in one interesting case in which one of the Lithium atoms captures it. 
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forming a LiH molecule. This theoretical formalism proves to be a powerful, 
effective and predictive tool for the analysis of non-adiabatic processes at the 
nanoscale. 
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1. Introduction 

The study of the interaction of charged particles with matter is a funda- 
mental area in modern physics, since these collisions are relevant for many 
fields of Science. Two relevant examples are radiation damage in biological 
tissues, and the stopping power of solids ^ essential in the design, for exam- 
ple, of fusion devices. In addition to the relatively old areas of ion-atom [ij] 
and ion-surface [21 collisions, the more recent intermediate discipline of ion- 
cluster collisions [3j has been developed in the last decades. 

In general, these processes - not only ion-cluster collisions, but all scatter- 
ing events of couples of nanoscaled objects, whether they are atoms, clusters 
or molecules, charged or not may trigger numerous processes, involving 
many of the degrees of freedom of the colliding projectiles: transfer of vibra- 
tional energy, transfer of electronic charge, ionization, fragmentation (some- 
times referred to as collision induced dissociation), recombination, fusion 
betweeen clusters, electronic excitations, etc. 

Perhaps the most important distinction that can be made is between adi- 
abatic and non adiabatic processes. The latter involve electronic excitations, 
and this fact implies, from a theoretical point of view, the necessity of a 
much more sophisticated method. If we constrain our consideration to ab 
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initio techniques, adiabatic processes could in principle be studied with stan- 
dard adiabatic first principles Molecular Dynamics (MD); however, whenever 
charge transfer, ionization, or simply electronic excitations play a role, some 
form of non-adiabatic treatment must be used. 

One such technique is Ehrenfest MD: it consists of assuming a classical 
approximation for the nuclei, that are coupled to the quantum electronic sys- 
tem. The latter is still a many-particle system whose out of equilibrium first- 
principles description is very demanding. It can be studied, however, ab ini- 
tio and non-adiabatically with the help of time-dependent density-functional 
theory (TDDFT) , that offers a good balance between computational 
eflFort and accuracy. This idea has received, among others, the names of 
Ehrenfest-TDDFT (E-TDDFT), non-adiabatic quantum MD (NA-QMD), or 
TDDFT-MD. 

The first application of the E-TDDFT equations was however in the realm 
of solid state physics [6], and in fact with the purpose of performing adia- 
batic MD in a diflFerent manner. Their use for non-adiabatic processes was 
pioneered by Saalmann and Schmidt jjl; a formal derivation of the model by 
Gross, Dobson and Petersilka can be found in Ref. 811. It has later been ap- 

n 

plied to various collision problems: ion-fuUerene collisions [9|, atom-sodium 
cluster collisions [10], charge transfer in atom-cluster collisions 



stopping power of protons or antiprotons in clusters Q or insulators 
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12|, the 



the excitation and ionization of molecules such as ethylene due to proton 
collisions [16], or the interaction of protons or heavier ions with carbon 
nanostructures or graphitic sheets Q. It may also be used to study 
laser-induced molecular or cluster dynamics in the high-field (but still not 
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relativistic) regime; some examples are Refs. flSl, l20|, l2l|, l22|, |23|, l2J, |25[ |. 

In this work, we have focused on the coUision of protons with the hthium 
tetramer. We investigate the feasibihty of using the E-TDDFT approach 
for identifiying the various possible reaction channels: these collisions are 
nothing else than chemical reactions with various possibles outcomes, that 
depend on the initial velocity, impact parameter, relative orientations, and 
even the initial vibrational state. Depending on the nature of the reactants 
and on their relative velocity, the reaction may be non-adiabatic, meaning 
that the electronic excited states play a role. A fully unconstrained first 
principles study of these chemical reactions in real time is far from possible ~ 
the complete study would require running over all possible initial velocities, 
orientations, etc - but the E-TDDFT model may provide some information 
for selected initial configurations. Moreover, the current experimental ad- 
vances in ultra-fast time-resolved observations of chemical reactions in real 
time demand parallel theoretical tools. 

In Section [2] we briefiy recall the theoretical framework and the compu- 
tational methodology; Section [3] describes the results; finally we summarize 
in Section [H 
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2. Methodology 

The Ehrefest MD scheme is defined by the following equations (atomic 
units are used hereafter): 

= H,{Rit)Mx,t) , (1) 
Mj^Rj{t) = -jdxip*{x,t)VjHeiRit)Mx,t) 

-V. E - ' (2) 

where (f{x^t) is the many- electron wavefunction, that depends on aU the 
electronic degrees of freedom, denoted x. It is governed by the electronic 
Hamiltonian HQ{R{t))^ which is determined by all the classical nuclear posi- 
tions R{t) = {Ri{t)^ . . . , RM{t)}' The motion of the nuclei is determined by 
the set of equations ([2]) - which are Newton's equations of motion for each 
nucleus J (characterized by a mass Mj and a charge Zj). The electronic 
Hamiltonian is given by: 

^ ^ 1 



j=i j<k 

-E - - • (3) 

The vector operators fj{j = l,...,N) are the electronic position opera- 
tors. This form of the Hamiltonian allows to write the force that acts on each 
nucleus solely in terms of the electronic density n(f,t), i.e., we can rewrite 
Eqs. [2] as: 

d^ - f 
Mj—Rj{t) = - d\n{r,t)VjvUr,R{t)), 

-V. E (4) 
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N 

1=1 

The electron-nucleus potential v^^ is given by: 



(f,t) = Mt)l5]5(f-f,)kW). (5) 



J \Rj{t)-r\ 

The possibility of computing the ionic forces solely in terms of the electronic 
density permits to use TDDFT, and propagate the proxy Kohn-Sham (KS) 
system of non-interacting electrons instead of the real one. Therefore, we no 
longer propagate Eq. [U but rather the time-dependent KS equations: 

z^(^,(f, t) = {-^ + ^KS [n] (f, t)}^^{r, t) = 1, . . . , N/2) (7) 

N/2 

n{f,t) = J2^\^iir,t)\'. (8) 

i=l 

Here, we have assumed an even number of electrons A^, and a spin-restricted 
configuration in which all the KS spatial orbitals (fi are doubly occupied. 
The non-interacting electrons move in the KS potential v^s [^] (^^ t) -> which is 
divided into the following terms: 



/nir' t) 
dV^^+^xeW(f,t). 



(9) 



The last term, v^c[^]{f^t)^ is the exchange and correlation potential, whose 
precise form is unknown and must be approximated. In this work we have 



chosen to use the simplest form, the adiabatic local-density approximation 
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We have used the octopus code to run these equations. The numerical 



details can be found in Refs 



271 . |28| ; here we will only summarize the essential 
aspects: In this TDDFT implementation, the wave functions, densities and 
potentials are discretized in a real space grid, instead of being expanded 



in basis sets. One important simplification arises from the use of pseudo- 
potentials (of the Troullier-Martins [29] type in this work), which substitute 
the nucleus and the core electrons of the system by a smooth set of local and 
non-local effective potentials, which are seen by the valence electrons. The 
Coulomb discontinuity disappears, and the valence orbit als no longer need to 
be orthonormal to the core ones. These facts allow for the use of a real space 
grid with a relatively large grid spacing (0.25A in this work). The molecule 
must then be placed in a simulation box, which is a sphere of 20A for the 
simulations that we present here. 

The first step in the simulation is the obtention of the ground-state of 
the Li4 cluster - both the lowest energy ionic geometry and the correspond- 
ing electronic structure. The former is a rhombic planar geometry (see 
Fig. [U the bond length is 2.92A). The proton is then placed 18A away 
from the centre of the cluster, and is given an initial relative velocity of 
0.02 a.u, which corresponds with a kinetic energy of 10 eV. The E-TDDFT 
equations are then propagated, and we have chosen a total simulation time 
T of 300 h/eV ^ 200 fs. In order to numerically propagate the equations, 
we use the velocity Verlet algorithm for time-stepping Newton's equations 



of motion, and an exponential midpoint rule for the TDKS equations j30|. 
Part of the electrons may leave the simulation box during the simulation - 
either because of ionization, or because some of the nuclei may also leave the 
simulation box, carrying away some electron density and this is accounted 
for by making use of absorbing boundaries. 
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Figure 2: Scattering trajectories for various (0, 6>) incidence angles. Atoms turn from red 
to white, and then to blue as time progresses. In this series, the d angle is kept fixed at 
7r/2. 
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Figure 3: Scattering trajectories for various (0, 0) incidence angles. Atoms turn from red 
to white, and then to blue as time progresses. In this series, the <p angle is kept fixed at 
7r/2. 
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3. Results 

An exhaustive study of this kind of colhsions would require the systematic 
variation of the foUowing initial conditions: the relative initial velocity of the 
colliding fragments, the vibrational state of the cluster, the incidence angles, 
and the impact parameter. Obviously, such a comprehensive study is out 
of the scope of any first principles methodology, and one must concentrate 
on a given parameter set and range. In this work, we concentrate on the 
incidence angles, since these are perhaps the variables whose small variations 
most easily produce different reaction outcomes. 

We have simulated the collision of the proton with the lithium tetramer 
at a relative velocity of 0.02 a.u, which corresponds with a kinetic energy of 
10 eV. This velocity was selected because it is sufficiently high to induce non- 
adiabatic effects (and indeed, we observed a small but non-negligible amount 
of ionization in most cases), but also not too high to induce too simple 
reactions: at higher velocities, we mainly observed "cluster transparency" 
(the proton passes through the cluster loosing part of its kinetic energy, but 
otherwise not altering its trajectory), followed in most cases by Coulomb 
explosion. 

The proton is directed towards the center of the cluster (the impact pa- 
rameter is therefore zero), and we have varied the incidence angles (0, 9) (see 
Fig. [1]), in order to study how the change in the incidence angle modifies the 
possible reaction outcome. We have grouped the simulations in two groups: 
in the first group, we have studied in-plane collisions, fixing at 7r/2 and 
varying the angle. These are shown in Fig. [21 In the second group, shown 
in Fig. [3l the cj) angle is fixed at = 7r/2, and we varied 9. Each collision 
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event is represented on a single picture; the complete trajectory is shown by 
superimposing all the snapshots (taken at intervals oi At ^ 5 fs), varying the 
color and tone of the atoms at different times to make the evolution evident. 

In the first series, it can be seen how in all cases except one, the proton 
is scattered away at varying angles. For one case, however (0 = ), one of 
the Lithium atoms captures it, and a LiH molecule is formed. Interestingly, 
the Lithium that captures the proton is the first one to collide with the 
proton, but the capture happens only after the proton is scattered away also 
from a second proton. This demonstrates an important fact: in order for the 
proton to be captured by the cluster or by any of the resulting fragments, 
its energy must be low. A large incidence energy, however, deos not rule 
out the possibility of a capture, because the proton may loose energy if 
it is successively scattered by more than one atom. Other than the LiH, 
two Li atoms associate and form a dimer. In other trajectories of the first 
series, the Lithium atoms reorganize in different manners: for example, it is 
interesting to see how the = collision results in Li2 plus two free Lithium 
atoms, whereas the apparently similar = 7r/2 case results in Lia plus only 
one Lithium atom. In the rest of the trajectories, the final configuration of 
positions (and velocities, which are not displayed) usually permits to predict 
what would be the final fragments. However, some cases are perhaps dubious, 
and longer simulation times would be needed. 

In the second series (Fig. [3j), the proton is scattered away in all cases. 
The reason is that the trajectories followed by the proton pass further away 
from any of the Lithium nuclei than in some of the in-plane cases. The first 
trajectory {9 = 0) is an example of cluster transparency: the trajectory of 
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the proton is unaltered, and the Lithium tetramer, although highly excited 
both vibrationally and electronically, remains intact. For ^ > 0, i.e., for non- 
orthogonal collisions, the proton leaves the cluster at a different angle. For 
> jq^^ the Li4 also retains its integrity, although strongly distorted. The 
other cases display different possible reaction channels. 

In all these plots we have not shown the electronic cloud. As an example, 
we display in Fig. [Hsix snapshots during the collision at (^f , f ) incidence 
angles. The analysis of the evolution of the charge may help to predict, at 
earlier times, whether or not two atoms are to remain associated: for example, 
the last snapshot (botton right) clearly shows the presence of charge between 
the two Lithium atoms on the left of the picture, and between the Lithium 
and the proton on top. This fact allows to infer with some confidence the 
formation of stable bonds (another possibility that would help in this task 
would be the use of the time-dependent electron localization function, which 
is also accessible with this methodology, see for example Ref. |3l|). 



4. Conclusions 

We have simulated the scattering of a proton with a lithium tetramer at 
moderate energies (10 eV), by making use of non-adibatic molecular dynam- 
ics based on TDDFT. An exhaustive analysis of this process would require a 
systematic variation of a large amount of initial conditions, such as the initial 
velocity of the colliding fragments, vibrational state of the target cluster or 
the impact parameter. This is unrealistic for any ab initio approach, but a 
careful selection of relevant parameters may provide interesting information: 
We have studied, for example, the different reaction channels that occur when 
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Figure 4: Six snapshot during the scattering trajectories for the (0, 6>) = (^f,f) in- 
cidence angles. The snapshots correspond with the times t = (top left), t = |T (top 
right), t = |T (middle left), t = |T (middle right), t = |T (botton left) and t 
(botton right). T ^ 200 fs. 



5 ^ 
6^ 
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the incidence angles of the proton change. Although long simulation times 
are needed in order to identify the final reaction products, these times are 
within the computational limits of the first principles approach that we have 
followed. In particular, we observed that for all the incidence angles consid- 
ered the proton is scattered away, except in one interesting case in which one 
of the Lithium atoms captures it, forming a LiH molecule. In any case, the 
resulting collisions reveal interesting information about the scattered proton 
and the different reorganization of the Li fragments after the collisions, pro- 
viding for each case self-explicative time-evolved trajectories. We conclude, 
therefore, that it is a suitable methodology to study the various reaction out- 
comes that result of proton-cluster collisions, or for any other non-adiabatic 
collisions in general. 
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